The landscape of two-proton radioactivity 
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Ground-state two-proton (2p) radioactivity is a decay mode found in isotopes of elements with 
even atomic numbers located beyond the two-proton drip line. So far, this exotic process has been 
experimentally observed in a few light and medium-mass nuclides with Z < 30. In this study, using 
state-of-the-art nuclear density functional theory, we globally analyze 2p radioactivity and for the 
first time identify 2p decay candidates in elements heavier than strontium. We predict a few cases 
where the competition between 2p emission and a decay may be observed. In nuclei above lead, the 
a decay mode is found to be dominating and no measurable candidates for the 2p radioactivity are 
expected. 

PACS numbers: 21.60.Jz, 23.50.-l-z, 23.60.+C, 21.10.Dr, 21.10.Tg 
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Introduction — With the impressive progress in map- 
ping new territories in the nuclear landscape, new phe- 
nomena emerge in rare isotopes with extreme proton-to- 
neutron imbalance. On the proton-rich side, due to the 
presence of the Coulomb barrier that has a confining ef- 
fect on the nuclconic density, relatively long-lived proton 
emitters exist beyond the proton drip line [l|-l3| . In recent 
decades, the phenomenon of proton emission from odd-Z 
nuclei developed into a powerful spectroscopic tool yield- 
ing a wealth of detailed structural information (sec the 
recent review [ij and references quoted therein) . In favor- 
able conditions unbound evcn-Z nuclei may undergo 
a simultaneous emission of two protons, i.e., exhibit 2p 
radioactivity. In such cases, due to proton pairing, the 
emission of a single proton is energetically forbidden or 
strongly suppressed. The ground-state 2p radioactivity 
was experimentally discovered in ^^Fe and, later on, 
in i^Mg 0], "^Ni 0], and ^^Zn 0. The interest in the 
phenomenon of 2p radioactivity has been boosted signifi- 
cantly by the measurement of proton-proton correlations 
m the decay of "^Fe ^ that has revealed the three-body 
character of the process and the sensitivity to the angu- 
lar momentum composition of the wave function. These 
findings were corroborated by the recent study of 2p cor- 
relations in the decay of ^Be resonances [llj. 

One may ask whether the ground-state 2p radioactiv- 
ity is limited to just a narrow range of light and medium 
mass nuclei or whether it can also be expected in heavy 
systems. No detailed predictions, however, have been 
made for elements heavier than strontium. Most of the 
previous theoretical estimates were focused on a rather 
narrow range of nuclei with 22 and 
aimed at identifying the best candidates for initial ex- 
perimental observations. Motivated by astrophysical ap- 
plications, these studies were subsequently extended to 
the region 30 < Z < 38 In almost all of these pa- 



pers, one- and two-proton separation energies of = —T 
nuclei were accurately determined (up to ~100keV) by 
calculating the Coulomb displacement energies in combi- 
nation with known experimental masses of mirror = T 
systems. The only exception is Ref. [l^ where self- 
consistent mean-field theory with various effective inter- 
actions was employed. 

The main objective of this work is to delineate for 
the first time the full landscape of 2p radioactivity. To 
this end, we use separation energies predicted by large- 
scale mass table calculations using state-of-the-art nu- 
clear density functional theory (DFT) [l3| with several 
Skyrme energy density functionals (EDFs) . In our global 
survey, we consider all even-Z elements with Z > 18. To 
estimate half-lives, we use two models of 2p emission: a 
direct-decay model and a diproton model. In addition, 
we take into account the competition between 2p emis- 
sion and a decay. Although our method is less precise 
than the approach based on Coulomb displacement ener- 
gies, it is well suitable for a large-scale, qualitative survey 
the 2p emission phenomenon undertaken in this study. 

Models — The nuclear binding energies B{Z, N) were 
obtained in the deformed DFT calculations of Refs. [isl . 
ITij using six effective Skyrme interaction models in the 
particle-hole channel (SkM* 0, SkP [2lLSLy4 SV- 
min [13, UNEDFO Q and UNEDFljli]) augmented 
by the density-dependent, zero range pairing term. The 
binding energies of even-even nuclei across the mass ta- 
ble were calculated by solving the self-consistent Hartree- 
Fock-Bogoliubov (HFB) equations using the solver HF- 
BTHO [26| . To approximately restore the particle num- 
ber symmetry broken in HFB, we used the variant of 
the Lipkin-Nogami scheme formulated in Ref. 27|. The 
binding energies for odd-iV isotopes were determined by 
adding computed average pairing gaps to the binding en- 
ergy of the corresponding zero-quasiparticle vacuum ob- 
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tained by averaging binding energies of even-even neigh- 
bors. Considering the uncertainties of current approaches 
to odd-even binding energy differences 128| , this is a rea- 
sonable procedure. 

The single-particle basis consisted of harmonic oscilla- 
tor states originating in 20 major oscillator shells. While 
the proton chemical potential Ap is positive for pro- 
ton unbound nuclei, the HFB results obtained with the 
discretized continuum are very stable in the considered 
range of binding energies. This is because the Coulomb 
barrier tends to confine the proton density in the nuclear 
interior and effectively pushes the continuum up in en- 
ergy [29I . [30j on the proton-rich side. As discussed in 
Ref. [18[, because of the Coulomb effect, the proton drip 
line lies relatively close to the valley of stability; hence, 
the associated model extrapolation error is small. Indeed 
all the models we use are very consistent when it comes 
to the prediction of the two-proton drip line, see Fig. [T] 

The half-lives for 2p decay were estimated using two 
simple models. The first, direct-decay model, results 
from the factorization of the decay amplitude into a prod- 
uct of two-body terms The removal of one proton 
leaves the core-|-p system in a state of energy Ep^ relative 
to the three-body decay threshold, and requires a trans- 
fer of orbital angular momentum Ip. The core-t-p sys- 
tem is taken here as the ground state of the one-proton 
daughter; hence, Ep = Q2p — Qp, where Q2p and Qp 
denote the decay energies for 2p and single-proton emis- 
sion, respectively. All calculations in this global survey 
were made with Ip = 0, i.e., assuming the fastest decay 
possible. In this way, we establish a limit of the least 
neutron-deficient nuclei decaying by the 2p emission. Wc 
note, however, that inclusion of larger values of angular 
momentum, in particular /p = 1, known to occur around 
Z = 28, would increase the number of predicted can- 
didates. The direct-decay widths were calculated using 
the version of the model given by Eq. (20) of Ref. [3j 
with the spectroscopic factor 9^ determined by compari- 
son with the experimentally established four 2p-emitters 
shown in Table [D Using the experimental separation en- 
ergies, the average value 0^ = 0.173 was obtained that 
gives a very reasonable agreement with experiment, see 
Table |IJ and has been used in subsequent calculations of 
half-lives: = ;iln2/r2p. 



TABLE I. Experimental partial 2p half-lives used to optimize 
the spectroscopic factors and the resulting predictions of the 
direct-decay and diproton models. In the direct model, Ip = Q 



was assumed. 


Nucleus 


Experiment 


direct 


diproton 


«Mg [7] 
^^Fe [10] 
"^Ni [8] 
''^Zn 


4.0(15) ps 
3.7(4) ms 
3.0l i J ms 
1.98io';Ii ms 


6.2ps 
1.1 ms 
6.8 ms 
1.0 ms 


12.3 ps 
8.7 ms 
5.3 ms 
0.8 ms 



The diproton model assumes that both protons leave 
the core nucleus as a correlated 2p pair with / = 0. 
Within this model 12 , l3| , the 2p-decay width is given by 
the Wentzel-Kramcrs-Brillouin expression. In our calcu- 
lations, the average diproton potential has been approx- 
imated by the 2V^(r), where Vp is the average proton 
potential containing the Woods-Saxon field in the Chep- 
urnov parametrization and the Coulomb term. (The 
results are fairly insensitive to the choice of the average 
potential [S]-) The diproton spectroscopic factor can 
be estimated in the cluster overlap approximation [l^ : 
= G'^[AI{A - 2)]2"e»2^ where = (2n)!/[22"(n!)2] 
|33| , is the proton overlap function, and n is the aver- 
age principal proton oscillator quantum number given by 
n « (3Z)i/3 _ 2 [3||. The value of = 0.015 was deter- 
mined by a optimization to the experimental half-lives 
of i^Mg, "Spe, ''^Ni, and ^"Zn. The values of half-lives for 
these nuclei predicted by the diproton model are given in 
Table H] they are consistent with the direct-decay model 
and the estimates of Refs. (31 . 



To determine the competition between 2p and a de- 
cay, the 2p decay half-lives were compared to a decay 
half-lives obtained from the global phenomenological ex- 
pression of Ref. . 

Selection criteria — The candidates for 2p decay were 
selected according to the energy criterion: 



!2p 



= -S2p > 



, Qp — Sp <c 0, 
B(Z, N) and S2p 



(1) 



where Sp = B{Z - l,iV) 
2, A^) — B{Z,N) ?a — 2Ap arc the one- and two-proton 
separation energies, respectively. The condition ([T]) cor- 
responds to true 2p decay as the simultaneous emission 
of two protons; the sequential emission of two protons is 
energetically impossible (see the inset in Fig. [T]). For the 
EDFs used in this work, the root-mean-square (RMS) de- 
viation from the experimental S2p values is typically less 
than 1 MeV. For instance, for UNEDFO and UNEDFl, it 
is 0.86 MeV and 0.79 MeV, respectively 

In addition to the energy constraint ([1]), we imposed 
the condition on 2p decay half-lives: 



lO^^s < T; 



2p 



< lO^^s, 



(2) 



which defines the feasibility of experimental observation 
of the 2p decay. The lower bound of 100 ns corresponds 
to the typical sensitivity limit of in-fiight, projectile- 
fragmentation techniques Q. The upper bound of 100 ms 
ensures that the 2p decay will not be dominated by (3 de- 
cay. (We note that the half-lives of the observed medium- 
mass 2p emitters are all in the range of several ms.) More- 
over, to eliminate the fast alpha emitters from our con- 
siderations, we only considered cases satisfying 



T2p < 10 -T^. 



(3) 



This condition guarantees that the 2p-decay branch is at 
least 10%. Of these candidates, to select the cases where 
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FIG. 1. (Color online) The landscape of ground-state 2p emitters. The mean two-proton drip line (thick black line) and 
its uncertainty (grey) were obtained in Ref. [3 by averaging the results of six interaction models. The known proton-rich 
even-even nuclei are marked by yellow squares, stable even-even nuclei by black squares, and known 2p emitters by stars. The 
current experimental reach for even-Z nuclei (including odd-A systems) [s^] is marked by a dotted line. The average lines 
Ni^vi^Z) of 2p emission for the diproton model (dashed line) and direct-decay model (dash-dotted line) are shown. The energetic 
condition ([1} for the true 2p decay is illustrated in the inset. 



the competition between 2p and a decay can be seen, we 
used the criterion 0.1 • T2p < < 10 • T2p, which ensures 
that the branching ratio for 2p or a decay is at least 10%. 

Results — For each model considered in this work, we 
selected candidates for 2p emission according to the im- 
posed criteria on lifetimes ^ and We define the 
model multiplicity ■m{Z,N) = fc if a nucleus {Z,N) is 
predicted by k models (/c = 1, . . . 6) to be a 2p emit- 
ter. The average path for the 2p emission in the {Z,N)- 
plane is given by Ni^^{Z), where - for a given element 
Z - the model-averaged neutron number is N^^iZ) = 
J2n N m{Z, ^)/J2n TTT-iZ, N), provided that at least one 
candidate has been found for this Z. Figure [T] shows the 
trajectories Nav{Z) for both the direct-decay and dipro- 
ton models. It is seen that (i) both ways of estimating 
2p half-lives give very similar predictions for the average 
path of 2p radioactivity and that (ii) this path quickly 
departs from the two-proton drip line with increasing 
atomic number. Furthermore, according to our calcu- 
lations, a decays wins over 2p emission above lead, so 
Z ^ 82 marks the upper border of the ground-state 2p 
emission landscape. The inter-model consistency for the 



predicted Q2p values along Ng^^f^Z) is good; namely, the 
RMS deviation for our six EDFs is typically 150 keV, i.e., 
well below the average deviation from experiment. 

Results of our survey are presented in more detail in 
Figure O We see that each element between nickel and 
lead has isotopes predicted to undergo 2p radioactivity. 
In the case of xenon {Z = 54), all 2p-decaying candidates 
are found to be dominated by a decay in the diproton 
model; -'^°^Xe is predicted to be a 2p emitter by the direct- 
decay model. For three light elements {Z = 20, 24, 26) no 
2p candidates are predicted because the calculated half- 
lives were shorter than the lower limit of condition ([2]) 
which is a consequence of our restriction to / = decay 
channel. We note that the observed 2p decay of ^^Fe is 
dominated by the I = 1 channel Q. While the nuclei 
54Zn, 59Ge, 63Se, 67Kr, and ^^Sr, discussed in are 
generally expected to meet the energy criterion ([1]) , their 
predicted Q2p values are too low to meet the lifetime 
criterion ([2]). In general, due to large uncertainties in the 
calculated half-lives due to uncertainties in Q2p [35[ , the 
estimated error on the predicted neutron number of a 2p 
emitter is A7V = 1. 
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FIG. 2. (Color online) The predictions of the direct-decay 
model (a) and diproton model (b) for the ground-state 2p ra- 
dioactivity. For each value of Z > 18, neutron numbers A'' 
of predicted proton emitters are shown relative to the aver- 
age two-proton drip line of Ref. [13 ] shown in Fig. [T] The 
model multiplicity m{Z, N) is indicated by the legend. The 
candidates for competing 2p and a decay are marked by stars. 
The current experimental reach of Fig. [1] is marked by a dot- 
ted line. 



In the region beyond ^^Zn, the predicted 2p-decay 
candidates which are closest to the current experimen- 
tal reach are ".ssQg (3^2), '^^^^^Se (2,1), ^eRr (3), and 
I02,i03fpg (3,2), where the numbers in brackets indicate 
the corresponding number of neutrons beyond the most 
neutron-deficient isotope known to date. All other cases 
are located by more than 3 neutrons from the present 
body of known isotopes. This distance is increasing with 
increasing atomic number and reaches 14 neutrons for 
i65pi3 which is predicted to be the 2p emitting lead iso- 
tope closest to the drip line. Other best candidates for 
ground-state 2p radioactivity in heavy nuclei, according 
to both direct-decay and diproton models, are: ^"^Zr, 
77Mo, siRu, 85pd, ii3ce, ii7Nd, i2iSm, ^^^^^^^Gd, "ODy, 

133-135gj. 138,139Yb 151,152Qg 154-156p^ g^-^^ 158,159jjg 



Two nuclei, -"^^^Pt and ^^^Hg, have been consistently pre- 
dicted to be 2p emitters in all models. 

In several cases a competition between 2p and a de- 
cay is predicted. The best candidates, marked in Fig. [2] 
by a star, are i^^Xe, los-noBa, i^'H^Cc, i27Gd, isipy, 

"4^145jjf^ aj^^ 147-149^^ rj,^^ j^^^jg- 144jjf and "8,149^ 

are predicted both in the direct-decay and diproton mod- 
els, but they are far from the line of the current exper- 
imental reach. The closest one is ^^-^Te, predicted in a 
diproton model with SkM* and SLy4 (Qsp « 3.3 MeV, 
Qa w 4.4 McV). (More recently optimized functionals 
SV-min and UNEDFl give Qzp « 2.75 MeV, 0„ w 
3.7 McV, i.e., much longer half-lives.) 

Conclusions — In this theoretical survey, based on the 
nuclear DFT, we quantified the landscape of ground-state 
2p radioactivity. To assess model-dependent extrapola- 
tions beyond the two-proton drip line, we applied six 
models based on Skyrme EDFs and two approaches to 
2p half-lives. Our results provide a consistent picture of 
the 2p radioactivity. Most importantly, we find that this 
decay mode is not an isolated phenomenon, limited to 
a narrow range of light and medium mass nuclei, but a 
typical feature for the proton-unbound isotopes with even 
atomic numbers. According to our calculations, almost 
all elements between argon and lead have 2p-decaying 
isotopes. The upper end of the 2p-decay territory is 
determined by a decay, which totally dominates above 
Z = 82. Unfortunately, most of the new candidates for 
the 2p radioactivity are located far beyond the current 
experimental reach. Only in two regions is the 2p decay 
mode predicted to occur close enough to be addressed 
by today's experiments. One ranges from germanium to 
krypton, and the second region is located just above tin. 
Other regions will have to wait for the facilities of the 
next generation. A confrontation of our predictions for 
heavier 2p emitters with the future data will be of great 
value for modeling of proton-unstable nuclei and improv- 
ing the nuclear EDF. 

Perhaps the most interesting are nuclei around ^°'^Te- 
^^"^Ba, in which the competition between 2p emission and 
a decay is predicted. The observation of these two de- 
cay modes in the same nucleus would provide an excel- 
lent test of nuclear structure models and a deeper un- 
derstanding of the dynamics of charged particle emission 
from nuclei. Finally, we note that all EDFs employed 
in our study, including the traditional ones (SkM* , SkP, 
SLy4) as well as the recently optimized ones (SV-min, 
UNEDFO, UNEDFl) yield a similar range of 2p radioac- 
tivity: while details for individual nuclei differ because of 
high sensitivity of 2p and a-decay half-lives to predicted 
Q- values, the global trends presented in this survey seem 
to be fairly robust. 
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